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1 Introduction 



When a system is governed by an autonomous nonlinear algebraic partial differential 
equation (PDE), it frequently admits some permanent profile structures such as 
fronts, pulses, kinks, etc [IB] , and usually these profiles are mathematically some 
single- valued solutions of the travelling wave reduction (x, t) — > x — ct of the PDE 
to an ordinary differential equation (ODE). The physical motivation of the present 
work is to find such solutions in closed form. Since this is a difficult mathematical 
problem, we restrict here to a simple case (a third order nonlinear ODE) and solve 
it completely. The method we used here is a refinement of Eremenko's method 
used in [3] as well as |5] and [6\ which is based on the local singularity analysis 
of the meromorphic solutions of the given differential equations as well as the zero 
distribution and growth rate of the meromorphic solutions by using Nevanlinna 
theory. This is a very powerful method. For example, it was used by Eremenko [5] to 
characterize all meromorphic traveling wave solutions of the Kuramoto-Sivashinsky 
(KS) equations. In fact, Eremenko showed that all the meromorphic traveling wave 
solutions of the KS equations belong to the class W (like Weierstrass) , which consists 
of elliptic functions and their successive degeneracies, i.e.: elliptic functions, rational 
functions of one exponential exp(kz), fceC and rational functions of z. 

In general, even if we know that the solutions belong to the class W, it is still 
difficult to find their explicit form. To overcome this problem, we shall apply the 
subequation method introduced in [13] and developed in [3]. In order to emphasize 
the method, we will choose a test equation according to the following criteria: 

1. to have a small differential order n, 

2. to have only nonrational Fuchs indices, apart from the ever present —1 index, 

3. to be of the form vS n ' = P with P a polynomial of its argu- 
ments, 

4. to have movable poles of order one, 

5. to be complete in the classical sense [15] (see details in [2, p. 122]) i.e. to 
include all admissible nondominant terms, 

The requirement for nonrational Fuchs indices sets n > 3. Let us take the 
complete autonomous third order polynomial ODE with simple poles, 

d u'" + diuu" + d 2 u' 2 + d^u 2 u' + d^u 4 

+C\u" + CiUv! + CsU 3 + C4V,' + C$U 2 + CqU + C7 = 0. (1) 
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This equation is indeed complete in the sense that it includes all polynomial terms 
having a singularity degree at most equal to four, as seen from the generating func- 
tion 

1 

(1 -tu){\ -iV)(l -tV) 

= 1 + ut + (u 2 + u')t 2 + (u 3 + uv! + u")t 3 

+ (u 4 + mV + u' 2 + uu")t 4 + 0(t 5 ). (2) 

Let us choose one particular set of dominant terms (the ones with coefficients dj, 
which have quadruple poles) so as to enforce from the beginning the condition that 
the Fuchs indices be nonrational. After setting C3 = without loss of generality, our 
test equation will be normalized as 

Cqu'" + 6u + ciu" + c 2 uu' + C4U 1 + c 5 u 2 + c 6 u + c 7 = 0, (3) 

Let u be a meromorphic solution of the ODE ([3]). We first check that if u has 
a movable pole at z = Zq, then u has only three distinct Laurent series expansions 
at Zq. Note that if z is a pole of u, it must be a simple pole. Therefore, in a 
neighbourhood of z = z , the Laurent series of the meromorphic solution u is of the 
form 

u(z) = U-i(z - ZqY 1 + U + Ui(z — Zq) H , M_i ^ 0. (4) 

Denote a any one of the cubic roots of Cq. Substituting the above Laurent series 
into the ODE ([3]) and balancing the leading terms, we obtain = a, and Uq = 
(—1c\a + c 2 a 2 )/(24c ). We are going to prove that there are at most three distinct 
Laurent series expansions at zq. If one linearizes the ODE (J3]) around the movable 
singularity z = zq [2J p. 114], the resulting linear ODE has the Fuchsian type at zq, 
and its three Fuchs indices r are defined by 

(r + l)(r 2 - 7r + 18) = 0. (5) 

Hence, the Fuchs indices are equal to r = —1, (7 ± \J — 23)/2. Because of the 
absence of any positive integer in the set of values of r, all other cofficients Ui are 
uniquely determined [2j p. 90] by the leading coefficient u~\. Hence, there are at 
most three meromorphic functions with poles at z = zq satisfying the ODE ()3]). 

We shall study the third order nonlinear differential equation ()3]) and show that 
all meromorphic solutions of this differential equation belong to the class W. More 
specifically, our results are the following. 

Theorem 1. If the ODE has a particular meromorphic solution u, then u belongs 
to the class W . Moreover, a necessary and sufficient condition for the ODE |3|) to 
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admit a particular meromorphic solution is to belong to the following list, 

c 2 

S 3a ■ ci, c 6 = arbitrary, c 2 = 0, c 5 = 0, c 7 = 0, c 4 = — — ; (6) 

IZCo 

c 2 

: c 5 , c 6 = arbitrary, c t = 0, c 2 = 0, c 4 = 0, c 7 = — — ; (7) 

I/O 

c 2 -12a 3 c 4 ci(c 2 + 36a 3 c 4 ) 

&2A '■ c 4 ,c 4 = arbitrary, c 2 = (J, C5 = — , c§ 



c 7 



4a 4 ' 144a 6 
^12a 3 c 4 -c 2 )(36a 3 c 4 -llc?)_ 



: Ci , c 2 = arbitrary, c 4 



1536a 8 

44c 2 + 8acic 2 — a 2 c 2 



144a 3 

32c 2 - 24ac 4 c 2 - 7a 2 c| {c x + ac 2 )(12c 2 + 6acic 2 + a 2 c^ 



c 5 = 77— a , °e 



48a 4 ' ° ~ 144a 6 

(4ci + 3ac 2 )(48c 2 + 20acic 2 + a 2 c\) 

° 7 ~ 55296a^ ; (9j 

S± : Ci, c 2 , c 4 , C5 = arbitrary, 

1152a 6 c 6 = -56c? + 60ac^c 2 - 18a 2 c x c\ + a 3 c^ + 288a 3 cic 4 

— 144a 4 c 2 c 4 — 96a 4 c 4 C5 + 48a 5 c 2 cs, 
2 13 3 2 a 8 c 7 = -176c 4 + 128ac 3 c 2 + 24a 2 c 2 c 2 - 32a 3 CiC 3 

+ 5a 4 c^ + 2688a 3 c 2 c 4 - 1536a 4 cic 2 c 4 + 96a 5 c^c 4 - 6912a 6 c^ 

+ 128a 4 c 2 c 5 - 512a 5 C!C 2 c 5 + 224a 6 c 2 c 5 

+ 4608a 7 c 4 c 5 + 2304a 8 c^. (10) 



We shall apply Eremenko's method [5] to prove the first part of Theorem 1. Here, 
we shall assume the readers are familiar with the standard terminology and results 
of Nevanlinna theory. The standard reference of this theory are and [121 E] ( see 
also [5] for a quick introduction). Our argument is slightly different from that of 
Eremenko and it makes use of the following version of Clunie's Lemma ([T21 Lemma 
2.4.2], see also [T7]). 

Lemma 1. Let f be a transcendental meromorphic solution of 

pP(zJ) = Q(zJ), 

where P(z, f) and Q(z, /) are polynomials in f and its derivatives with meromorphic 
coefficients {a\\X G 1} such that m(r, a\) = S(r, /) for all A G I. If the total degree 
of Q(z, f) as a polynomial in f and its derivatives is less than or equal to n, then 

m{r,P(r, /)) = S(r, f). 
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Now let u be a function meromorphic in the complex plane which satisfies the 
above ODE ([3]). If u is rational, then we are done. So suppose u is a transcendental 
meromorphic solution of equation ([3]), then we have 

— 6u = CqU 1 " + Ciu" + C2UII + C4I1' + C5M 2 + CqU + C7, (11) 

Take / = u,P — u,n — 3 and apply Clunie's lemma (Lemma 1) to the above 
equation, we conclude that m(r, u) = S(r, u) and hence (1 — o(l))T(r, u) = N(r, u). 
We claim that u must have infinitely many poles. Assume it is not the case, then 
N(r,u) = O(logr). Therefore, T(r,u) = O(logr) which is impossible as u is tran- 
scendental. 

Secondly, we prove that if u is a transcendental meromorphic solution, then u 
is a periodic function. Recall that there are at most three meromorphic functions 
with poles at z = z satisfying the ODE (J3j). Now let Zj,j = 1, 2, 3, • • • be the poles 
of u(z), then the functions Wj(z) = u(z + Zj — Zq) are meromorphic solutions of the 
ODE fl3]) with a pole at z . Thus some of them must be equal. Consequently, u is 
a periodic function. 

Without loss generality, we may assume that u has a period of 27ri. Let D = 
{z : < Imz < 2tt}. If u has more than three poles in D, then by the previous 
argument, we can conclude that u is periodic in D and hence it is indeed an elliptic 
function and we are done. 

Now suppose u has at most three poles in D. Since u is a periodic function with 
period 27ri, we have N(r,u) = O(r), as r — > 00. It follows from (1 — o(l))T(r, u) = 
N(r,u) that T(r,u) = 0(r). By Nevanlinna's First Fundamental Theorem, we 
know that for any a G C, N(r, 1/(m — a)) = 0(r) as r — > 00. By the periodicity of 
u, we conclude that u take each a finitely many times in D. Hence, the function 
R(z) = u(\nz) is a single- valued analytic function in the punctured plane C\{0} and 
takes each a G C finitely many times. It follows that is a removable singularity 
of R and R can then be extended to a meromorphic function on C. Hence R is a 
rational function as it takes each complex number finitely many times. Therefore, 
u(z) = R(e z ) belongs to the class W and this completes the proof of the first part 
of Theorem 1. 

Remark. From the above proof, we notice that if u is an elliptic solution, then u 
has at most three (simple) poles in each fundamental polygon Q. Recall that the 
residue of u at any pole must be one of a, wo, co 2 a where u is the cubic root of unity. 
Since the sum of the residues of all the poles in any fundamental polygon Q is zero, 
u must have three distint simple poles in Q and hence we have three distinct Laurent 
series at z - 

Remark. If u(z) = R(e kz ) where R is some rational function, then R has at most 
three (simple) poles in C\{0}. We are going to show that R cannot have a pole at 
0. Suppose we write u(z) = R{Z) = r /Z n + Y^h=i r H ~ %i) + q(Z), where q is 
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a polynomial in Z = e kz . Substituting u(z) = R(Z) into ODE ([3]) and letting Z 
tend to infinity, we can conclude that q equals to some constant C . Now letting Z 
tend to 0, we can deduce that r = 0. Hence, u(z) = T, i= ir i /(e kz -Z i )+C, where 
Zi, C e C. Finally, if u is rational, then u will have at most three (simple) poles in 
C. Similarly, we can show that u must be of the form Yli=i r i/i z ~~ z %) + Cj where 
r„C6C. ' 

2 Explicit solutions in the class W 

Let us determine the constraints on the coefficients Cj of fl3]) for meromorphic solu- 
tions to exist, and let us determine all these meromorphic solutions in closed form. 
According to section [H these solutions are necessarily elliptic or degenerate of ellip- 
tic (i.e. rational in one exponential e kz , k G C or rational in z), i.e. they belong to 
the class W. 

If the meromorphic solution is elliptic, by a classical theorem, the sum of the 
residues of the three Laurent series for u, Eq. 03]), must vanish, and similarly for any 
rational function of u, u', u" . These necessary conditions pl[| are first established in 
section 12.11 

If the solution is elliptic, one knows the elliptic orders of u and u', they are 
respectively equal to three (three simple poles) and six (three double poles). There- 
fore, by a classical theorem of Briot and Bouquet [H p. 277], [TJ part II, chap. IX 
p. 329], [9, p. 424] the elliptic solution obeys a first order algebraic equation whose 
degree in u' is the order of u (three) and degree in u is the order of u' (six), 

m 2m— 2k 

F(u, v!) = ^ ^2 a jyk u J u' k = 0, a 0jm ^ 0, (12) 

k=0 j=0 

with m = 3. The complex constants a^, with ao, m ^ 0, are then determined by the 
algorithm presented in [13], i.e. by requiring each of the three Laurent series (]4]) to 
obey ( TT2|) . The search for all third degree subequations (|T2|) obeyed by the three 
Laurent series is performed in section 12.21 

As to those solutions of (EJ) which are degenerate of elliptic, they also obey a first 
order equation (TT2|) . whose degree m is at most three. Because of the singularity 
structure of ()3]) (three distinct Laurent series), any m-th degree subequation, 1 < 
m < 3, must have m distinct Laurent series. The search for all second or first degree 
subequations f|T2l is performed in sections (12. 3p and (12.41) . 

Let us first establish all these first order subequations. Their general solution 
may be either singlevalued (and hence in class W) or multivalued. The explicit inte- 
gration of the singlevalued subset will provide as a final output all the meromorphic 
solutions of ([3]) in closed form. 
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2.1 Residue conditions 

If ([3]) admits an elliptic solution, it is necessary that, for any rational function of u 
and its derivatives, the sum of the residues inside a period parallelogram be zero, 



VA: e N Vn e N res = °- 



(13) 



i=i 



The first conditions are 

k = 0,n = 2: c 2 = 

k = 0, n = 3 
fc = 0, n = 5 



u 



0. 71 

1, 7i 



C 4 



(14) 



12a 3 ' 
c 4 c 5 = 0, 

C1C7 = 0, 

(c 6 (c^ - 128c 7 ) = if ci = 0), (c 7 (c 3 + 36a^c 6 ) = if c x 0). 



When the computation is limited to < 4, n < 10, this defines three and only 
three distinct sets of fixed coefficients for a possible elliptic solution, 



c 2 = 0, ci = 0, c 4 = 0, c 6 7^ 0, c 7 = 
c 2 = 0, ci = 0, c 4 = 0, c 6 = 0, 



c 2 = 0, ci 7^ 0, c 4 



12a 3 ' 



c 5 = 0, c 7 = 0. 



(15) 
(16) 

(17) 



2.2 Subequations of degree three 

Denoting Uk, k = 1,2, 3, cubic roots of unity, each such subequation has the neces- 
sary form 

F(u,u') = —{uiau' + u 2 )(u! 2 au' + u 2 )(u 3 au + u 2 ) 



+ biu'^u + 6 2 u'u 3 + 6 3 u 5 + 6 4 u' 2 + b 5 u'u 2 + b§u 
+ 6 7 m'u + b$u 3 + 69?/ + 6f,u 2 + b a u + b c + b 

with all ujj distinct and the additional condition to be irreducible. 



0. 



(18) 



The first order third degree subequation is precisely defined as 

F(u, u') = -a 3 u' 3 - u 6 + btu^u + b 2 u'u 3 + 6 3 m 5 + b 4 u' 2 + b 5 u'u 2 + b 6 u* 

+ b-ju'u + b 8 u 3 + b 9 u + b b u 2 + b a u + b c + b = 0. (19) 



The algorithm [13] to compute the coefficients bk is to substitute u by one of the 
Laurent series (jl]), which makes the right hand side of fll9p become a Laurent series 



+00 



F(u,u') 



TFAz-zoy- 6 , 



/ j j 



(20) 



then to solve the infinite set of equations 

Va Vj : Fj = 0. (21) 

The practical resolution is as follows. First, the 21 equations Fj = 0,j = 0,...,6 
define a linear system for the which admits a unique solution and generates six 
nonlinear constraints among the six c&. By considering slightly more equations in 
(l2~Pj) (in the present case, going to j = 8 is enough), the set of nonlinear constraints 
among the c^'s admits exactly two solutions, and all the remaining equations Fj = 
identically vanish. These two solutions are 



r 2 

S 3a '■ Ci, Cq = arbitrary, c 2 = 0, c 5 = 0, c 7 = 0, c 4 



12a 

2 

4 



(22) 
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S-sb '■ C5, C6 = arbitrary, c\ = 0, c 2 = 0, C4 = 0, C7 

and they are identical to the two residue conditions (fT7j) and f TTS]) . 
The corresponding subequations have genus one 

(au + Ak lU ) 2 {au' - 2km) + (u 3 + 20kf + k 6 ) 2 = 0, c x = 12a 2 h, c 6 = 4k 6 , (23) 
(au'f + (m 3 - 3k 2 5 u + k 6 ) 2 = 0, c 5 = c 6 = 4fc 6 . (24) 

The method to integrate them [Tj §249 p. 393] is to build a birational transformation 
to the canonical equation of Weierstrass 

p' 2 = A(p - ei)(p - e 2 )(p - e 3 ) = 4p 3 - - 9s- (25) 

To do that, it proves convenient to introduce one of the roots e of the cubic poly- 
nomial of u(x) appearing as a square in (123"]) and fT24]) . i.e. to redefine k 6 by the 
respective relations 

e 3 + 20k 3 + k 6 = and e 3 - 3k 2 e + k 6 = 0. (26) 

The subequation ( 12~4|) is one of the five first order binomial equations of Briot 
and Bouquet [2, p. 122], its general solution is classical 

1 _ p'{z-z ,g 2 ,g 3 )-A (e 2 - k 2 f(e 2 - 4k 2 ) 



92 = 0, g 3 



u-e JVj ' ^ ' *° 243a 6 

2(eg - k 2 ) 2 e (e 2 - k 2 ) 

Nl ~ 3a3 ' ^ " 3a3 ' (27) 

The subequation ff23l has been integrated by Briot and Bouquet [l, §250 p. 395] 
by introducing a function w defined by 

X£ 3 — 6 3 

aw + 4k x u = Q-w, (28) 

u-e 
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then by establishing the birational tranformation 



w 



au' + Ak\U 



u z + e u + 



2 • 



U 



—3aww' — e w 3 + 6kiw 2 + 2e 
2(w 3 + 1) : 



finally by integrating the ODE for w, 

2ki A 4Jfei(Jfe?-eg) eg - 20egfc? - 8fcf 

W — 1 — , 92 — 7T~A > 93 — 



17a 6 

<,:; 27 9i '" . .4 ■ — /i - 



e p — B ' 3a 4 

(8A; 3 + e 3 ) 3 e 3 ^ _ _e|+8fc| 

3a 2 ' 



27a 
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a" 



(29) 



(30) 



More generally, birational transformations from (u,u') to (p,p') are obtained 
with an algorithm due to Poincare, implemented for instance by the command 
Weierstrassf orm of the computer algebra package algcurves [TT] . 



2.3 Subequations of degree two 

Let us define the second degree subequation as 

F(u, u') = a 2 u' 2 — au 2 u + u 4 + b^u'u + b^u 3 + b^u' + b 2 u 2 + b\U + bo = 0, (31) 

with the additional condition to be irreducible. Computations similar to those 
mentioned in section 12.21 provide two solutions, 



S2A '■ Ci,C4 = arbitrary, C2 = 0, C5 
Ci(c\ + 36a 3 c 4 ) 



12a 3 c 4 



4a 4 



c 6 
C7 



144a 6 

;i2a 3 c 4 -c 2 )(36a 3 c 4 -llc 2 ) 



1536a 8 

1 9 7 9 ?9 

u = v — — , ci = —3a fci, c 4 = 2ao + -0%, 



(32) 



2 L2\ 2 

it — tr N 



av — 



+ -(v + 6)(v - b)(v - hf = 0, 6^0, 
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and 



S*2_b : Ci,c 2 = arbitrary, c 4 



44c 2 + 8acic 2 — a l c^ 



-32cf - 24acic 2 - 7a 2 cj 
48a 4 

(ci + ac 2 ) (12c 2 + 6aciC 2 + a c: 

144a 6 ' 
(4ci + 3ac 2 )(48cf+ 20acic 2 + a 2 c|) 



144a 3 



c 5 



c 6 



C7 



(33) 



55296a 7 

u = v + - + — -g, c 2 = -2— + 6a6, 
4 12al a 



2\ 2 



at> — 



+ ^{v + bf{v -b) =0, 6^0. 



For k 2 ^ b 2 , the point transformation 



1 

u = ki H , w 

w 



1 



fci + b ki — b 
maps the ODE ([32]) to the Riccati ODE 

b 2 



+ N A 



iV 2 



(A; 2 - 6 2 )^ 



(34) 



aNX' - MX - 



A(k 2 - b 2 ) 



(A 2 + 1) = 0, M 



3b 2 



A(k 2 -b 2 )' 



(35) 



whose general solution is a Mobius function of one exponential so that v is a rational 
function of one exponential. 

For k\ = b 2 , i.e. for instance for k\ = —b, the ODE fl33|) integrates as 



2b 



6 + — , w = 1 + 3 (l + e 



b(z-z )/(2a) 



W 



(36) 



2.4 Subequations of degree one 

These first degree subequations 

F(u, u') = aw' + u 2 + 6iM + 6q = 0, 



(37) 
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are determined by requiring their vanishing when u is the Laurent series (J4]). This 
results in 

' Si : c%, c 2 , c 4 , C5 = arbitrary, 
2ci — ac 2 

61 " 12a 2 ' 

_ 44cf - 32acic 2 + 5a 2 c\ - 144a 3 c 4 + 144a 4 c 5 

= 1152a 4 ' 

1152a 6 c 6 = -56c? + 60ac^c 2 - I8a 2 cic 2 2 + a 3 c 3 2 + 288a 3 cic 4 

— 144a 4 c 2 C4 — 96a 4 ciC5 + 48a 5 c 2 cs, 
2 13 3 2 a 8 c 7 = -176c 4 + 128ac?c 2 + 2Aa 2 c\cl - 32a 3 cic 3 , 

+ 5a 4 c^ + 2688a 3 c 2 c 4 - 1536a 4 CiC 2 c 4 + 96a 5 cf,c 4 

- 6912a 6 c 2 + 128a 4 c 2 c 5 - 512a 5 ciC 2 c 5 + 224a 6 c 2 c 5 
+ 4608a 7 c 4 c 5 + 2304a 8 c 2 . 

The solution of this Riccati equation is either a rational function of one expo- 
nential or a rational function, 



(3* 



u 



bi k k. . 
-j + a-coth-(z- z ) 

bi a 



2a 2 



(39) 



2 ^ z — z 



, bf - 46 = 0. 
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